Hydrodynamic Interactions in Two Dimensions 
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We measure hydrodynamic interactions between colloidal particles confined in a thin sheet of fluid. 
The reduced dimensionality, compared to a bulk fluid, increases dramatically the range of couplings. 
Using optical tweezers we force a two body system along the eigenmodes of the mobility tensor 
and find that eigen-mobilities change logarithmically with particle separation. At a hundred radii 
distance, the mobilities for correlated and ant i- correlated motions differ by a factor of two, whereas in 
bulk fluids, they would be practically indistinguishable. We derive the two dimensional counterpart 
of the Oseen hydrodynamic tensor which quantitatively reproduces the observed behavior. These 
results highlight the importance of dimensionality for transport and interactions in colloidal systems 
and proteins in biological membranes. 



When governing a world of lower dimensionality, the 
laws of physics give rise to intriguing phenomena. The 
reduced number of spatial dimensions usually results in 
stronger and longer ranged correlations. One example 
of this is long range Coulomb correlations and electron- 
lattice interactions which give rise to peculiar electronic 
and structural phase transitions in system of low dimen- 
sionality [H. In a similar way, fluid flow propagators, 
mediating hydrodynamic interactions between suspended 
bodies, have logarithmic tails in two dimensions [2[, giv- 
ing rise to strong dynamical correlations 

When two particles suspended in a viscous fluid ap- 
proach each other, propagation of fluid flow results in hy- 
drodynamic interactions. As a consequence, a particle's 
motion resulting from a given stimuli, strongly depend 
on the entire instantaneous spatial configuration. Hydro- 
dynamic interactions tend to favour correlated motions 
where every particle can move within the "slip-stream" 
of its neighbours [4] . This has consequences on many col- 
lective phenomena such as colloidal aggregation and gel 
formation where attractive interactions push to modify 
interparticle distances [5[. Moreover particles are always 
subject to stochastic thermal forces whose effects may 
be strongly influenced by hydrodynamic interactions, as 
happens in polymer dynamics 0, or for protein confor- 
mational changes 0, S] • In three dimensional bulk flu- 
ids, the strength of hydrodynamic coupling decays as the 
inverse interparticle distance. This long range charac- 
ter makes hydrodynamic interactions quite effective in 
determining dynamical behavior and poses a number of 
numerical and theoretical challenges to the physical mod- 
eling of such phenomena. 

The situation is even more dramatic when fluid flow is 
restricted to two dimensions, such as in cell membranes 
or soap films, and flow fields decay logarithmically with 
distance. It is important to note that spatially confin- 



ing the particles to two dimensions is not sufficient to 
observe these long range tails. It is essential that the 
momentum flow is similarly restricted on a two dimen- 
sional plane. To this aim it is crucial that the bounding 
fluid has a much smaller viscosity than the film itself, a 
solid boundary would lead instead to hydrodynamic in- 
teractions decaying faster than in the 3D case [13] • 

The problem of two dimensional, single particle mobil- 
ities has received a lot of attention due to its relevance 
for Brownian motion in cellular membranes However, 
much less is known about the equally non trivial and im- 
portant role of many body hydrodynamic interactions in 
two dimensional systems. Using video microscopy, Che- 
ung and coworkers Q observed long ranged spatial corre- 
lations between spheres floating in a free standing liquid 
film. Manybody effects were found to be significant even 
at very small concentrations, practically precluding the 
possibility of isolating hydrodynamic effects from indirect 
bulk measurements. 

In this Letter we demonstrate how optical microma- 
nipulation allows the direct observation of hydrodynamic 
interactions between an isolated pair of colloidal particles 
suspended in a free standing liquid film. Optical tweez- 
ers allow the positioning of two beads at varying separa- 
tions, isolated from boundaries and other particles. The 
tweezers can also be used to simultaneously drive parti- 
cles along any direction and directly probe the full mobil- 
ity tensor. As a result, we observe a dramatic increase in 
the interaction range of hydrodynamic interactions due 
to the reduced dimensionality. On the other hand, we 
derive a two dimensional expression for the Oseen hy- 
drodynamic tensor. Once the relevant far field boundary 
conditions are recognized, we can quantitatively describe 
experimental observations. 

We are interested in determining particle motions un- 
der the action of an external force field. To the lowest 
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order in hydro dynamic interaction (large enough inter- 
particle distances), it is sufficient to know single particle 
quantities and in particular the mobility b and the far 
field Green's tensor G(r) describing flow propagation: 



= bFi 
u(r) = G(r - R- 
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where u(r) is the flow field produced by a force applied 
to a particle located in R^. 

In a many body system the velocity of the i-th particle 
will be (within the superposition approximation) the sum 
of two contributions: the speed that it would have in the 
absence of other particles plus an "ambient" velocity ob- 
tained as the sum of all the fluid velocities independently 
produced at R^ by other particles located at Rj. 



R.; = bFi 
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In bulk three dimensional fluids, b is the inverse Stokes 
drag (l/Gnfia) while G is the Oseen tensor [S.] which, in 
cartesian coordinates, reads: 



(4) 



If we neglect stresses produced by the fluid bound- 
ing the film (air in our case) then the variations of flow 
properties across the film are negligible and dynamics is 
governed by a two dimensional Stokes equation. In other 
words, we only allow for in plane momentum flow and ne- 
glect contributions from the film surfaces. Working with 
relatively slow flows at the micron scale, the Reynolds 
number will be negligibly small allowing to express the 
momentum balance equation in the form of a two dimen- 
sional Stokes equation: 



/iV^u(r) - Vp(r) = -Y5{v)/h 



(5) 



When complemented by the incompressibility condi- 
tion V • u = 0, equation ([5]) can be rewritten as a bihar- 
monic equation and solved for the propagator: 
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(6) 



where we have neglected terms, whose coefficients 
depend on particle size and shape and boundary condi- 
tions on its surface. Such terms soon become negligible 
when moving away from the particle and play no role 
for hydrodynamic interactions between sufficiently sepa- 
rated particles. However, these terms determine the ac- 
tual value of velocity on the particle boundary and there- 
fore the single particle mobility b. The logarithmic term 



appearing in (|6]) precludes the possibility of determin- 
ing the integration constant L by imposing a vanishing 
velocity at infinity. The divergence of flow field signals 
the presence of a length scale beyond which some of the 
assumed approximations fails. 

There are three approximations involved in the deriva- 
tion of Eq. ([5]): 1) infinite size of the film, 2) negligible 
inertia, 3) negligible viscous drag on the interfaces. For 
each of them there is a length scale beyond which the 
solution in (|6|) is not self consistent with the assumed 
approximation. For the first approximation, the length 
scale is clearly the actual size of the film Li. Following 
P we can impose a frictionless boundary on a ring of 
radius Li and obtain the expression (|6]) for the propaga- 
tor where L = Li. For a stick boundary condition one 
gets L = L\l ^Je [IjJ. In the second case, the flow field 
propagated by (|6|) fails to satisfy the negligible inertia ap- 
proximation when the distance r is of order L2 = l-i/ pU^ 
where p is the fluid density and U is the typical particle 
speed. However, inertial terms in the Stokes equation 
can be a ppr oximately taken into account by the Oseen 
method [l2[ and obtain a short distance (r <C L2) ex- 
pression for the propagator which is given by (|6]) with 
L = 4exp[l — 7]L2 — 6.IL2 (7 is the Euler-Mascheroni 
constant). At last, one can obtain from Saffman so- 
lution [9] that the momentum flow through film inter- 
faces cannot be neglected for distances of the order of 
L3 = hp/p' where h is the film thickness and /i, p' are 
respectively the viscosities of the film and that of the 
bounding fluid. However, even in that case, expression 
(|6]) remains valid in a neighborood r <^ Ls when we re- 
place L with exp[l/2 — j]Ls — O.9L3 . When dealing 
with mesoscopic systems however, due to the very low 
Reynolds number involved, L2 is usually very large, a 
few hundred meters in our case. L3 is also quite big 
when the bounding fluid has a much lower viscosity than 
the film fluid, as in the present case where L3 ~ 0.1m. 
Therefore, we expect that hydrodynamic interactions in 
free standing liquid films are usually dominated by the 
finite size of the membrane. The following discussion will 
remain valid even in other situations, provided one uses 
the relevant length scale L. 

We can arrange all R^ in a single 2N dimensional vec- 
tor and introduce the 2N x 2N hydrodynamic mobility 
tensor H so that (|3]) reads: 



H(R) F 
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Hf{K) = S.jS^^b + (1 - S,j)G^^{K, 



(7) 

R.) (8) 



We now have a two dimensional expression for the 
propagator ([6|) that can be used in (|8]) to get the many 
body mobility tensor as a function of the two parameters 
b and L. We will check the validity of our expression 
for H by observing the dynamical behavior of two col- 
loidal particles confined in a free-standing liquid film. In 
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FIG. 1: (Color online) Trapping geometry. Two beads of 
radius a are optically trapped at a distance r in a liquid film 
of tfiickness h. 



a two body system, if we choose the x axis along the 
joining line, x and y dynamics are naturally decoupled 
for symmetry reasons. We can then decompose H into 
two 2x2 mobility tensors operating on the subspaces of 
X and y coordinates. Any 2x2 symmetrical matrix is 
diagonal in the coordinate system of the two eigenvectors 
(1,1) and (1,-1), corresponding respectively to a rigid 
translation and a stretching motion. Without losing any 
generality we can then assume that our mobility tensor 
is fully characterized by its four eigenvalues. Using ex- 
pression in ([8]) we can diagonalize the mobility tensor 
in the two body case and obtain the four eigenvalues as 
a function of interparticle distance r: 
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FIG. 2: Relaxation of eigen-coordinates. For each of the 
selected distances, we report the time evolution of the four 
eigen-coordinates relaxing to equilibrium after the applied 
perturbation. Solid lines are exponential fits. Dashed line 
is the average, single particle dynamics. 

As in the three dimensional case, hydrodynamic inter- 



actions produce a splitting in the spectrum of mobilities. 
The splitting is symmetric about the average single parti- 
cle mobility and is larger for parallel than for perpendicu- 
lar forcing. In three dimensions, the entity of the splitting 
decays as the inverse separation falling below 10% when 
particles distance grows beyond ten radii. The depen- 
dence on distance, in the two dimensional case, occurs 
only through a logarithmic term, which makes hydrody- 
namic interactions practically unavoidable. 

In order to directly observe 2D hydrodynamic inter- 
actions, latex beads (2 fim diameter) are dispersed in 
a water-glycerol mixture with 0.2% wt SDS surfactant 
added. A thin film is obtained by sweeping the solution 
on a square frame (6 mm side) of nylon wires (60 /im 
thickness) [13]. Glycerol increases viscosity and slows 
down both drainage and evaporation, resulting in longer 
lived films. Starting with a 50% wt water /glycerol film 
and then heating to evaporate most of the water, we can 
obtain a very viscous film with a few micron thickness. 
We measured the film thickness h at the beginning and 
at the end of the reported experiment to be 3.9 /im. Par- 
ticles are imaged by a 40x NA 0.75 objective of an in- 
verted optical microscope (Nikon TE2000-U). The same 
objective is used to focus the laser beam (A=532 nm) 
diffracted off a spatial light modulator (SLM, Holoeye 
LCR-250Q) into two, dynamically reconfigurable, optical 
traps 15 1 . Axial confinement is by capillary force on 
the top surface of the liquid film (Fig. [1]). 

In contrast to previous studies of hydrodynamic in- 
teractions in 3D 

i, H [13. where eigen-mobilities 



where extracted from the correlated fluctuations of opti- 
cally trapped particles, we choose to measure the eigen- 
mobilities by directly exciting the four eigenmodes. For 
each of the four eigenmodes, we calculate two holograms 
producing two sets of traps slightly displaced along the 
eigencoordinate. The two holograms are alternatively 
displayed onto the SLM. The SLM is fully refreshed with 
the new hologram on a time scale (50 ms) which is much 
faster than the time scale of particle dynamics. The two 
particles will be then displaced from the new equilibrium 
positions, along the selected eigencoordinate. In our ex- 
periment, external forces are of optical and stochastic 



origin: 



-kSKi + Si 



(11) 



with k the trap strength and SHi = — the ith 
particle displacement from trap center R^. Assuming 
small displacements compared to interparticle distances 
and averaging over Brownian motion, we get from ([71 [TT]): 



((5R) = -m(R°) • (^R) 



(12) 



If the initial configuration corresponds to a small dis- 
placement along the nth eigenmode ((5R(0)) = ee^, then 
(fT2|) has the solution: 
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FIG. 3: Eigen-mobilities. The four eigen-mobilities of a two 
particle system arranged at different particles separations. To 
correct for variations in trapping power, for each distance 
mobilities have been normalized to their average. 



{SH) = eex.-p[—kXnt]en 



(13) 



and we can directly obtain the corresponding mobility 
An by monitoring the amplitude of the mode relaxing 
to equilibrium. We choose to normalize the eigenvec- 
tors such that the corresponding eigencoordinates give 
the center of mass position, for the rigid modes, and the 
half distance (along x or y) for stretching modes. 

We choose ten logarithmically spaced interparticle sep- 
arations between 5 jam and 100 /im and drive the two 
particles back and forth (eight times) along each of the 
four eigenmodes. Particle coordinates where digitally ex- 
tracted from video frames at 144 Hz. Eigencoordinates 
where then computed and averaged over the eight itera- 
tions. Fig. [2] shows the time evolution of the four eigen- 
coordinates at four selected interparticle distances. The 
two X and y rigid motions are much more mobile than the 
corresponding stretching motions. This behavior remains 
clearly visible up the highest investigated distance (100 
particle radii). The four eigen-mobilities corresponding 
to the four probed modes can be extracted by fitting 
the eigen-coordinate dynamics to the exponential law in 
(p!3|) . To correct for variations in trap strength /c, for each 
distance, we normalize the four obtained decay rates to 
their average value kb. The relative eigen mobilities are 
shown in Fig. [3] as a function of the particle separation. 
The strength of hydrodynamic coupling, reflected in the 
splitting of mobilities, decays logarithmically slow with 
distance. Still at a separation of 100 radii , particles move 
twice as fast when forced along the same direction rather 
than in the opposite. At the same large separation, three 
dimensional mobilities would only differ by 1%. 

The four data sets in Fig. [3] can be very well fitted by 
(|9|) leaving L and the adimensional mobility 6* = Anfihb 



as the only free parameters. We obtain as best fit pa- 
rameters 6* =8.7 and L = 2.1 mm, the corresponding 
fitting curves are shown as straight lines. As expected, 
the relevant length scale L is determined by the film fi- 
nite size. In particular a sticky boundary condition on 
a ring inscribed in the film frame (Li = 3 mm) would 
give an L = 1.8 mm. For the same boundary condi- 
tion, the single particle mobility b can be calculated for 
a cylinder of height /i, the film thickness [11]. Using 
our particle radius a as the cylinder radius we obtain 
6* = log(Li/a) — 1 = 7.0 which compares reasonably well 
with the corresponding fitted value. We do not expect 
the two values to be in better agreement since the mobil- 
ity b depends on the details of the boundary conditions 
on the particle surface. We anticipate that, when moving 
to multiparticle systems, the mobility of long wavelength 
eigenmodes will diverge linearly with the number of par- 
ticles, rather than logarithmically as in the 3D case. As a 
first consequence of that, we could transport, at the same 
speed, any number of beads, using the same total amount 
of laser power. Moreover, the crossover to underdamped 
propagating modes on a linear chain of trapped particles 
predicted in [17|, could be experimentally verified much 
more easily. 

We have directly measured hydrodynamic interactions 
between colloidal particles in a sheet of viscous fluid. The 
reduced dimensionality, compared to the bulk 3D fluids, 
results in stronger and longer ranged hydrodynamic cou- 
plings which are quantitatively reproduced by a two di- 
mensional version of Oseen hydrodynamic tensor. The 
observed interactions constitute a general model for diffu- 
sion and interactions of proteins in biological membranes. 
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